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The algorithm is a derandomization of the ITERATEDRADON algorithm of Clarkson et al. 
(International Journal of Computational Geometry and Applications 6 (3) (1996) 357-377) 
and is guaranteed to terminate with an 2(1/d7)-center. Moreover, it returns a polynomial- 
time checkable proof of the approximation guarantee, despite the coNP-completeness of 
testing centerpoints in general. We also explore the use of higher order Tverberg partitions 
to improve the running time of the deterministic algorithm and improve the approximation 
guarantee for the randomized algorithm. In particular, we show how to improve the 
Q(1/d?)-center of the ITERATEDRADON algorithm to 2 /d=) for a cost of O((rd)*) in 
time for any integer r > 1. 
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1. Introduction 


A centerpoint of a set S of n points in R@ is a point c such that every closed half-space containing c also contains at 


least TAI points of S. Intuitively, every hyperplane through c divides S into roughly equal parts. The Centerpoint Theorem 


asserts that every set of n points in R@ has a centerpoint. This was first established by a theorem of Rado [14] for general 
measures and by Danzer et al. [4] for the special case of point sets. 

A centerpoint is a natural generalization of the median to higher dimensions. In fact, a centerpoint is an approximate 
median of the set f(S) for every linear projection f : S —> R. They are used as robust estimators in statistics, because they 
are invariant under affine transformation and robust to outliers [5]. They are also used in mesh partitioning [7]. 

The Tukey depth of a point x is the minimum number of points of S lying in a half-space that contains x. Thus, a 
centerpoint is a point with Tukey depth at least TT More generally, c is a -center if it has depth at least £n. We consider 


the problem of finding an approximate centerpoint, giving a deterministic algorithm for computing a zay center. 


The existence of centerpoints can be proven directly from either of two classic theorems of convexity theory, Helly’s 
Theorem and Tverberg’s Theorem. In Section 3, we discuss how these two proofs of the centerpoint theorem lead to different 
perspectives for designing algorithms for computing approximate centerpoints. 

The exact complexity of computing centerpoints in higher dimensions is not known. The decision problem of testing 
if a given point is a centerpoint is coNP-complete [16]. However, a simple corollary of Tverberg’s Theorem guarantees the 
existence of a subset of centerpoints, call them Tverberg points, that admit polynomial-time checkable proofs. Moreover, 
testing if a point is a Tverberg point is NP-complete [16]. In this case, the decision problem is well understood but sheds 
little light on the hardness of the search problem of actually finding a centerpoint. 
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The fastest known algorithm for computing a centerpoint of S C R@ is due to Chan [1] and computes a -center in 
time O(n‘~’) in expectation, where £ is the maximum achievable for the set S. Such a 8-center is also known as a Tukey 
median. 

The ITERATEDRADON algorithm of Clarkson et al. was the first algorithm that computes an approximate centerpoint in 
time sub-exponential in d [2]. The algorithm computes an 2(1/d*)-center with high probability. Section 3 describes how 
this algorithm resembles the proof of the Centerpoint Theorem via Helly’s Theorem. 

The main operation in the ITERATEDRADON algorithm is to replace sets of points by their Radon point, a point in the 
common intersection of the convex hull of two disjoint subsets. Radon’s Theorem guarantees the existence of such a point. 
Tverberg’s Theorem is a generalization of Radon’s Theorem that guarantees a common intersection for a larger collection of 
subsets. 

In this paper, we use the intuition from Tverberg’s Theorem to iteratively construct a proof of depth for an approximate 
centerpoint. The result is a new approximation algorithm, ITERATEDTVERBERG, that derandomizes the ITERATEDRADON algo- 
rithm of Clarkson et al. In Section 4, we prove that the ITERATEDTVERBERG algorithm produces an 2(1/d*)-center in time 
sub-exponential in d. Moreover, the center comes with a polynomial-time checkable proof of its depth. 

We elaborate on this intuition in Section 5, showing how solving larger sub-problems can be used to speed up the 
running time of the deterministic algorithm and to improve the approximation ratio of the randomized version. 


2. Related work 


Centerpoints are the most well known definition of a geometric median [6]. Like many such medians, it can be computed 
via linear programming and the problem of finding a “best” centerpoint can be written as a maximum feasible subsystem 
problem (see [5] for a survey of computational aspects of data depth). The linear constraints are just the set of hyperplanes 
intersecting d+ 1 affinely independent points of S. There are O(n‘+') of these. Consequently, any linear programming 
method will require time n°, limiting their usefulness to low-dimensional instances. 

In the plane, centerpoints can be computed in linear time [9]. Several algorithms are known to compute centerpoints 
in R? in O(n? polylogn) time [3,13]. The best known algorithm for d > 3 is due to Chan and runs in O(n4—') expected 
time [1]. Chan’s algorithm computes the deepest possible centerpoint, also known as a Tukey median. He conjectures that 
the O(n*—') running time is optimal for this problem in the algebraic decision tree model. However, the exact complexity 
of computing centerpoints is not known. In particular, it is not known if it is possible to compute a centerpoint in time 
polynomial in n and d. 

Several approximation algorithms for centerpoints exist in the literature. Many approaches using random sampling are 
known [11,16,2]. Verbarg showed that for dense point sets, the mean is a good approximate centerpoint [18]; it is a B-center 
where £ depends on the density. 

The only previously known algorithm to compute an approximate centerpoint in time sub-exponential in d is the ITERAT- 
EDRADON algorithm of Clarkson et al. [2]. The ITERATEDRADON algorithm returns an 2(1/d7)-center with high probability in 
time polynomial in n and d. ITERATEDRADON is a Monte Carlo algorithm and there is no way known to verify that the point 
returned by the algorithm is indeed an 2(1/d*)-center. The inner loop of ITERATEDRADON depends on the following classic 
theorem [15]. 


Theorem 2.1 (Radon’s Theorem, 1921). Given n > d + 1 points S C Rf, there exists a partition (U, U) of S such that conv(U) N 
conv(U) Æ Ø. 


We call a partition of d+ 2 points as described in the theorem, a Radon partition, and we call a point in the intersection, 
a Radon point. 

The simplest version of the ITERATEDRADON algorithm works as follows. Build a balanced (d + 2)-ary tree of height h. 
Fill in the leaves with points from the input set S by sampling them uniformly at random. Each interior node of the tree 
is filled in with the Radon point of its children. Clarkson et al. show that a height of h = Ign is needed to compute an 
92(1/d?)-center with high probability, resulting in a running time that is O((dlgn)!84) [2]. Thus it is sub-exponential in d 
but not polynomial. 

Several different modifications to the basic ITERATEDRADON algorithm are presented in [2]. In one version, the running 
time is reduced to O(poly(n,d)) by reusing the Radon points the same way the input points are reused, but this bound 
hides the dependence on 4, the probability guarantee that the output is correct. Other variations of the algorithm incorpo- 
rate linear programming to solve larger subproblems. These variants are able to produce {$ -centers with some constant 
probability. 

It is not known how to derandomize the other versions of ITERATEDRADON. We do, however, explore a variation on the 
linear programming based solution. We present a new way to leverage these larger subproblems and analyze the impact on 
both the randomized and the deterministic algorithms (see Section 5). 
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Fig. 1. Radon points for R? and R?. Tverberg points for r =3 in R? and R°. 


3. Two proofs of the Centerpoint Theorem 


Theorem 3.1 (The Centerpoint Theorem (Rado [14], Danzer et al. [4])). Given a set of n points S c R, there exists a centerpoint c € Rf 
such that every closed half-space containing c also contains at least Ta points of S. 


The Centerpoint Theorem is most often presented as an easy consequence of Helly’s Theorem. It is also possible to 
prove the existence of centerpoints via Tverberg’s Theorem (see [12] for a comprehensive treatment). The relationship 
between these two proofs gives insight into the relationship between the ITERATEDRADON algorithm and its derandomization 
presented in this paper. 


Theorem 3.2 (Helly’s Theorem [8], 1913). Given a collection of compact, convex sets X1,..., Xn C Ri. If every d +1 of these sets have 
a common intersection, then the whole collection has a common intersection. 


The Centerpoint Theorem follows from Helly’s Theorem as follows. Consider the set H of all open half-spaces that contain 


more than -4 points of S. For each such half-space h € H, let Pa denote conv(S Mh), a compact, convex set containing the 


same points of S as h. Any d + 1 of the half-spaces have a common intersection at one of the points of S, and thus every 
d+ 1 of the Py’s also have a common intersection. We apply Helly’s Theorem to the family of sets Py = {Pp | h € H}. The 
common intersection guaranteed by Helly’s Theorem is exactly the set of all centerpoints. 

The most common elementary proof of Helly’s Theorem makes extensive use of Radon’s Theorem, despite that Helly’s 
Theorem technically came first (though published second). The proof first considers the case where there are only d + 2 
sets. The hypothesis of the theorem implies the existence of d + 2 points, each taken from the common intersection of 
d-+1 of the sets. The Radon point of these d+ 2 points satisfies the conclusion of the theorem. The proof for n > d+ 2 
sets uses induction. The inductive step again considers a set of points taken from each of the common intersections of 
n—1 sets, and shows the Radon point of this set satisfies the theorem. Unraveling this induction into an algorithm, we 
arrive at something very much like the ITERATEDRADON algorithm except with far too many Radon point computations to 
be computationally feasible. This is why ITERATEDRADON uses random sampling used to avoid paying for the combinatorial 
blowup in the number of sets. The ITERATEDTVERBERG algorithm we present shows how to make this sampling deterministic 
(Section 4). 

The Centerpoint Theorem can also be proven via Tverberg’s generalization of Radon’s Theorem. 


Theorem 3.3 (Tverberg’s Theorem [17], 1966). Given (d + 1)(r — 1) + 1 points S c Rf, there exists a partition of S into S1,..., Sy, 
such that ()j_, conv(S;) 4 Ø. 


Observe that Radon’s Theorem is a special case of Tverberg’s Theorem when r = 2. 

Say that a point c is a Tverberg point if it is in the common intersection of the convex hulls in the Tverberg partition. 
Then, setting r= Ed yields a Tverberg point contained in the convex hull of [za] pairwise disjoint subsets of S. Any 
half-space containing c must also contain at least one point from each of the subsets and therefore, c is a centerpoint. 

Observe that the centerpoints guaranteed by Tverberg’s Theorem come equipped with a polynomial-time checkable proof. 
Given the partition, we need only verify that the point is in the convex hull of each part. If any part in the partition has 
more than d + 1 points then by Carathéodory’s Theorem, there is a subset of size d+ 1 that contains the Tverberg point is 
its convex hull. We may therefore assume the convex hulls are simplices of dimension at most d, so checking can be done 
quickly. The key insight in derandomizing the ITERATEDRADON algorithm is to actively construct these Tverberg partitions for 
the intermediate points used in the algorithm. 


4. Derandomizing the ITERATEDRADON algorithm 
4.1. Short proofs of depth 


Given a point p € Rf, a proof that p has Tukey depth (at least) r with respect to S is a collection of r disjoint subsets 
U1,..., Ur C S such that p € N= conv(U;). If there is such a proof, we say p has depth r. For any input point s € S, the 
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set {s} is a proof of depth 1. A point and its proof of depth is denoted by (point, proof). We say that a collection of proofs 
are disjoint if their sets are pairwise disjoint. 

The ITERATEDTIVERBERG algorithm is very similar to the ITERATEDRADON algorithm. The key difference is that each succes- 
sive approximation carries with it a proof of its depth. When we combine d + 2 points of depth r into a Radon point c, we 
can rearrange the proofs to get a new proof that c has depth 2r as shown in the following lemma. 


Lemma 4.1. Given a set P of d + 2 points with disjoint proofs of depth r, the Radon point of P has depth at least 2r. 


Proof. Let (P1, P2) be the Radon partition for P, and let c be the Radon point. For each p; € P, order the parts in the proof 
partition of p; and call the jth part U; j. We build a proof that c has depth at least 2r with parts of the form Ler U; j for 
ke {1,2} and j € {1,...,r}. These sets are disjoint by our hypothesis. Since each p; € Pg is contained in conv(U;, j), it follows 
that Pk C conv (U piep, Ui j) and also, by standard convexity arguments, conv(Pk) © conv(U pier, Ui, j). By the definition of 
Radon points, c € conv(P;). Therefore, c € conv(P,) © conv(Up,<p, U; j) for all 2r choices of k € {1,2} and j €{1,...,r}. 


The preceding lemma implies a simple deterministic algorithm for computing an approximate centerpoint. Construct a 
(d+ 2)-ary tree with n leaves. Store the points of S in the leaves. In each interior node of the tree, store the Radon point of 
its children. The height of the tree is logy, .n, so Lemma 4.1 implies that the depth of the root is 2!8+2" = O(n!/'84+2))_ Not 
too shabby for such a simple algorithm, but the depth guarantee of the output is only sub-linear in n. To get an 2(1/d?)- 
center, we need to find a way to build this tree higher, and in order to do that, we need more leaves. The following lemma 
gives a hint as to where we can look for more points to put in the leaves. 


Lemma 4.2. If there is a proof that a point p has depth r, there exists such a proof that contains at most r(d + 1) points of S. 


Proof. Let P;,...,P; be the sets in the proof for p. By Carathéodory’s Theorem, there exists a subset P; C P; of at most 
d+ 1 points such that p € conv(P}). So, the sets P{,..., P; is the desired proof of the correct size. 


We refer to this economizing of proofs as pruning. A standard proof of Carathéodory’s Theorem implies a reasonable 
algorithm for pruning a proof. Starting with the point of interest c as a convex combination of k points (where k > d+ 1), 
find a new convex combination that uses one less point by subtracting a multiple of an affine dependence of the points. 
The choice of the affine dependences does not matter so we may find one by solving a linear system on any subset of d+ 2 
points in O(d?) time. We then iterate this procedure until we have c as a convex combination of just d + 1 points. 

As in Lemma 4.1, a set to be pruned is the union of one set from each proof of each point in one part of a Radon 
partition. If the proofs to be combined were already pruned, and thus contained only sets of size d+ 1, then the total 
number of points in the combined sets is at most (d + 1)*. Consequently, the pruned set can be found by after O (d?) 
iterations. The entire pruning operation takes O (d°) time. 


4.2. The deterministic algorithm 


We are now ready to describe the ITERATEDTVERBERG algorithm. We make buckets By to store (point, proof) pairs for 
points of depth 2°. At the start of the algorithm, the input points s € S are placed in bucket Bọ with the trivial proof of 
depth 1. The main loop of the algorithm takes any bucket Be with a least d+ 2 points qi,...,qq42 and combines them as in 
Lemma 4.1 to form a proof of depth 2°+! for the Radon point of q1,...,qg42. We then prune the proof as in Lemma 4.2 and 
push it into bucket Be,,. The loop terminates when we get a point in bucket Bz, where z = Era A single iteration 


of the main loop is illustrated in Fig. 2 and pseudocode is provided in Algorithm 1. 
4.3. Theoretical guarantees 


Theorem 4.3. The ITERATEDTVERBERG algorithm returns a point of depth at least lza l 


Proof. Clearly, if the algorithm terminates then it returns a point of the desired depth and its corresponding proof. Suppose 
for contradiction that the algorithm does not terminate. Then, at some time t, it is impossible to pick i because every bucket 
has fewer than d + 2 points. Let h be the maximum such that Bp is nonempty at time t. For any point p in a bucket B;, the 
proof for p contains at most 2'(d+ 1) input points. Since all of the proofs in all of the buckets are disjoint, we have 


h 
X zd+1? >n. (1) 
i=0 
Therefore, 
n 
gh+1 — (2) 


7 (d41)2’ 
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a) p > b) 
c) d) 


e) 


Fig. 2. The ITERATEDTVERBERG algorithm: (a) Sets of d +2 points are divided into Radon partitions. (b) d+ 2 Radon points are combined into a second-order 
Radon partition. (c) A proof polygon is formed by taking the convex hull of two subsets, one from each of the Radon points in the second order partition. 
(d) The proof polygon is pruned to a simplex. (e) All of the proof polygons before the pruning phase. (f) The proof simplices after the pruning. 


f) \ 


Algorithm 1 ITERATEDTVERBERG 
Input: $< R4:|S|=n 
Initialize empty stacks Bo,..., Bigi 


Push (s, {s}) to Bo for each se S 
while B,,;__» | is empty do 


its =a 
2(d+1)2 


2(d+1)2 
Initialize proof to be an empty stack 


Let £ be the max such that By_, has at least d+ 2 points 
Pop d+2 points qy,..., qa+2 from Be—1 
Let (c, {U1, U2} = RADON(q1,..., da+2) 
for k=1,2 do 
for i=1,...,2°-! do 
Let S;; be the ith part of the proof for qj 
Let X= Ugjeur Sij 
Let X’ =PRUNE(c, X) 
Push X’ to proof 
Push ((s, {s})) to Bo for each xe X \ X’ 
end for 
end for 
Push (c, proof) to Be 
end while 
return (c, proof) from Bigr n} 
2(d+1)2 


and thus, 
h n 


FeLi (3) 


However, this is impossible because it implies that h > lg ay and therefore Bigp_» _j is nonempty. This would have 


2(d+1)2 
caused the algorithm to break out of the while loop and terminate, contrary to our assumption. 


652 G.L. Miller, D.R. Sheehy / Computational Geometry 43 (2010) 647-654 


Theorem 4.4. The ITERATEDTVERBERG algorithm requires n° “8% time. 


Proof. The main loop of the algorithm increases the size of some By by one and decreases the size of Bg_; by d+2. It 
requires a Radon partition computation and 2° pruning operations that dominate the running time of the inner loop. 

As mentioned previously, the Radon partition can be found in O(d*) time, and the pruning takes O(d°) time. Since 

n 
L<lg MENA these two operations require O (d°) + 2 zarn? 0 (d?) = O (nd?) time. 

We will count the number of iterations of the main loop, charging the iterations for a point that was popped from a 
bucket to the point created in that iteration. Let the cost of point in bucket Bẹ be defined inductively as one more than the 
sum of the costs of the d + 2 points that have been popped from Be_;. Since the cost of a point in Bẹ is determined by £, 
it can be written as a simple recurrence, 


te =(d+2)te1 +1. (4) 
In the base case to = 0, because input points have trivial proofs of depth 1. It follows that 

te < (d+2)°. (5) 
Let Tg be the total number of iterations of the main loop given buckets of sizes B = {|Bo|,...,|Bm|}. The value of Tg 


can be computed exactly by summing the costs of the points remaining in buckets. 


m 
Ts =) |Bilti. 
i=1 
At termination, no bucket has size greater then d+ 1. So, if 6 is the set of final bucket sizes, then 


m 
Tg < X @4+ ti < tmy- (6) 
i=1 
Eqs. (5) and (6) imply the Tg < (d + 2)"*!. By Theorem 4.3, the final depth is arD” so we are interested in m = 


[lg xan The total number of iterations of the main loop can now be bounded as follows. 


Tg < (d+ 2) zam! (7) 
7 ( i i (8) 
(d+ 1)2 
= nO (lsd) (9) 
Since each iteration of the main loop requires O (nd?) time and the number of iterations is n°“°S®, the total running 
time is n? 084), 


5. Leveraging larger subproblems 


The Radon point of a set of d+ 2 points has Tukey depth at least 2 and, thus, is a centerpoint for that set. Both the 
ITERATEDRADON algorithm and the ITERATEDTVERBERG algorithm replace sets of d+ 2 points with their centerpoint using 
Radon’s Theorem. In this section we address the effect on these algorithms if we instead solve larger subproblems. Rather 
than combining points in sets of d+ 2, we look at sets of size (d+ 1)(r — 1) + 1 for some fixed r. The centerpoint of such 
a set has Tukey depth at least r. It is not known how to solve these larger problems in time sub-exponential in d. However, 
if n is large and d is not too large, it may be feasible to solve subproblems in O (df) time even though O (nô) is prohibitive. 
In fact, we will see that this can even lead to an asymptotic speedup for the ITERATEDTVERBERG algorithm. 


5.1. Improving the approximation for the ITERATEDRADON algorithm 


The Radon point of d+ 2 points is a centerpoint of the subset. Consider the following modified version of the ITERATED- 
RADON algorithm. 

We can run the same iterative algorithm as before except using r-partitions instead of 2-partitions. In fact, it is not 
necessary to keep around the partition, it actually suffices just to find any centerpoint of (d+ 1)(r — 1) + 1 points at 
each round. We can go through the analysis from [2] and see the impact of r in the depth guarantee of the approximate 
centerpoint. 

The analysis works by looking at any projection of the point set to the line. We compute the probability f(x) that the 
tree of iterations with height h returns a center of depth at most nx. Without loss of generality, the projections of the points 


of S land on a £,..., 1. It follows that fo(x) < x. The center of the center will be x where x is such that f(x) is very small. 
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At each iteration, the centerpoint is at least r deep in the projection. There are (+1) CLER) choices for the r points less 


than the centerpoint in the projection. The probability of the union over all these choices is bounded from above by the 
sum of the individual probabilities, and thus, 


ai 1(r—1)+ ‘) 
r 


fr x) < fn—1(x)’. (10) 


-1 
Let B = Coe) : 


r 


Fa < BO fas)" 11) 
< BTH (BT fr) 12) 
< BBB fo)” 13) 
<p 14) 
rh 
<er ( x ) 15) 
Bri 


Now, since 6 = O((rd)~"), we can choose x smaller than O((rd)=1) and the probability f,(x) vanishes as desired. So, 
even for a choice of r = 3, we can improve the depth guarantee of the resulting centerpoint by an O(/d) factor. 


5.2. Speeding up the ITERATEDTVERBERG algorithm 


The same trick of solving larger subproblems can speed up the running time of the deterministic algorithm. Tverberg’s 
Theorem guarantees the existence of a partition of S into r sets whose convex hulls have a common intersection as long 
as |S| > (d+ 1)(r — 1) + 1. Say T(r) is the time required to compute a Tverberg partition into r parts. To the best of our 
knowledge, nothing better than brute force is known for computing Tverberg partitions for r > 2. 

A slight modification to the ITERATEDTVERBERG algorithm to use Tverberg r-partitions instead of Radon partitions results 
in an n'8"/T(r) speedup. Thus, for n large enough, we get an asymptotic speedup. 

The modified algorithm simply makes recursive calls on sets of [n/r] points and combines them in sets of (r — 1)(d + 
1) +1. The analysis is virtually identical to the original version except we give up a factor of r/2 in the depth of the output. 
As for the running time, the new algorithm now has a recursion tree with higher fan out and the resulting running time is 
O((d + 2) "T(r)) = O(n'8E+2)/I8rT(r)), 


6. Conclusions and open problems 


We have presented the ITERATEDIVERBERG algorithm, the first algorithm that deterministically computes an approximate 
centerpoint in time sub-exponential in d. By combining intuition from both Helly’s Theorem and Tverberg’s Theorem, our 
method sheds an interesting new light on the problem of computing approximate centerpoints. It still remains open whether 
it is possible to compute approximate centerpoints deterministically in time polynomial in n and d. We conjecture that it is. 

We also extended both our algorithm and the ITERATEDRADON algorithm by looking at the impact of solving larger 
subproblems. One consequence of this work is that any new results on quickly computing centerpoints for small point sets 
can be used to improve these algorithms. Currently, it is not known how to compute centerpoints of more than 2d + 2 
points in time polynomial in d. However, we conjecture that computing the centerpoint of 2d + 3 points in R? is NP-hard. 

In the ITERATEDRADON algorithm, it was not clear if the factor of d appearing in the output depth approximation ratio 
was intrinsic to the problem or merely an artifact of the analysis. This same factor of d shows up in our analysis of the 
ITERATEDTVERBERG algorithm in a completely different way, implying that perhaps it is intrinsic to the problem. Finding the 
optimal approximation ratio that can be achieved by a deterministic algorithm in sub-exponential time remains an open 
problem. 

It is also open as to whether the ITERATEDTVERBERG algorithm can be significantly sped up by doing the pruning step 
lazily. In such a variant, the “extra” points that get pruned out would be points with proofs of their own, possibly allowing 
the reuse of more computation. Thus far we have been unable to analyze the lazy pruning version of the algorithm. 

The reader is directed to Kalai’s survey for many other interesting problems related to Tverberg points [10]. 

The computation of centerpoints draws a compelling correspondence between fundamental theorems in convexity theory, 
Helly’s Theorem and Tverberg’s Theorem, and fundamental complexity classes of NP and coNP. It is our hope that future 
work will further elucidate this correspondence. 
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